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The interaction of an edge dislocation made of half the superconducting 



> 

o 

^1 ' plane with a magnetic interlayer vortex is considered within the framework 

Q\ , of the Lawrence-Doniach model with negative as well as positive Josephson 



interlayer coupling. In the first case the binding energy of the vortex and the 



dislocation has been calculated by employing a variational procedure. The 

O 

current distribution around the bound vortex turns out to be asymmetric. In 
the second case the dislocation carries a spontaneous magnetic half-vortex, 
whose binding energy with the dislocation turns out to be infinite. The half- 
vortex energy has been calculated by the same variational procedure. Impli- 
cations of the possible presence of such half- vortices for the properties of high 
temperature superconductors are discussed. 
PACS numbers: 74.72.-h, 74.20.De, 74.20.Hi 



1 Introduction 

Many properties of high temperature superconducting (HTS) materials are well un- 
derstood within the framework of Lawrence-Doniach model [1,2] assuming Josephson 
type coupling between the Cu-0 superconducting (SC) layers. Recently in [3] it has 
been proven experimentally that the current-voltage and microwave characteristics 
of these compounds exhibit features specific to a stack of Josephson junctions (JJ), 
each being 10-i-15 A thick in the c-direction [3]. To date, the LD model was widely 
employed (see [4-7]) to describe the behavior of interlayer magnetic vortices in an 
ideal layered structure. 

In [8] Annett has suggested that the Josephson interlayer critical current J c might 
be negative, if the BCS pairing mechanism is confined within each layer and does 
not work between the layers. In this case the ground SC state of the LD model 
is characterized by the phase shift ir of the order parameter in adjacent layers. To 
distinguish such a state from the conventional one [1,2] and following the terminology 
of [9], the former and the latter will be named below as 7r-superconducting ground 
state (-7T-SC) and O-superconducting ground state (0-SC), respectively. It is worth 
noting that several mechanisms may account for the presence of 7T-SC (or J c < 0). 
As was suggested in [9], if the tunneling between two metals is controlled by spin-flip 
processes, the Josephson critical current between these metals should be negative. 
The other mechanism causing J c < [10], relies on the possibility that the tunneling 
could occur indirectly through some intermediate state where Coulomb repulsion 
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excludes a double occupancy of such a state. We will not further discuss here the 
extent to which these mechanisms are pertinent to the interlayer tunneling in HTS 
compounds. Instead, we will consider how the existence of 7r-SC, if any, may be 
recognized in specific features of vortex dynamics. 

The 0-SC and 7T-SC appear to be physically indistinguishable in an ideal infinite 
crystal unless a screw or edge dislocation is present [11]. Then in the case of the n- 
SC these should carry a spontaneous half-vortex (HV) as a consequence of the phase 
mismatch ir, while traveling around the dislocation line. Therefore an investigation 
of interaction between vortices and dislocations can distinguish these two SC ground 
states. 

In what follows we will consider the interaction of the interlayer magnetic vortex 
(which has no normal core) with the edge linear dislocation made of half a SC plane 
in the two cases (0-SC and 7r-SC) of the LD model. 

2 Model 

We will employ the LD model [1,2]. As was discussed in [2,6], because of weak 
interlayer coupling, the constant amplitude approximation for the in-plane order 
parameters can be adopted in order to describe interlayer vortices. Thus, the LD free 
energy functional becomes dependent on the phases <£>„ and the in-layer components 
of the vector potential A n , n labels the layers. The edge dislocation can be modeled 
as a line z = 0, x = so that for x > the conducting half plane n = is missing 
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(FIG.l). To incorporate it into the LD model it is useful to separate the contribution 
of the n = plane. Consequently, the free energy functional without the dislocation 
can be written as [2] 

/s + OO —OO 1 

dh \j2(K n + U+) + £ (K n + U~) + — dz (curl A) 2 

L n=l n=-l ° n 

+ K + U } , (1) 
where the notations are introduced as follows: 



K n = _ (^o) 2 — (V$ n + — A„) 2 , 

Kb <fo 

Un = 7^M l ~ ^C0S($„ - $ n±1 -Xn,n±l)] , 

C/q = ^J c [2-z/cos($ -$-i -X0-1) -^cos($i - $ - Xi,o)] , 



27r rn+l 

Xn,n+1 = / 

ipO Jn 



Here, the superconducting layers are infinitesimally thin and separated by insulating 
layers of thickness s; the coordinates r = (x, y) lie in the planes, and z is directed 
perpendicular to the layers; ip stands for the unit flux, X ab is the in-plane London 
penetration length, J c is the critical interlayer current; V = (d x ,d y ),A = (A x ,A y ). 
One can see that only the last two terms of (1) contain $ . 

The parameter v is introduced here to distinguish between the two possible states 
of the LD model: 0-SC [y = 1) [2] and vr-SC ( v = -1) [8]. In the format of (1) 
these two states can be converted into each other by the phase shift $ n — > $ n + irn. 
However, the presence of a dislocation makes such a conversion impossible. 

To describe a dislocation in the n = layer (FIG.l) one should replace the 



Josephson (potential) energy Uq by the following expression: 

n I 2-i/cos($o-$-i-Xo,-i)-*'cos($i-$o-Xi,o), x < 
Uod = — J c { 

[ 1 -v cos -x_i,i), a:>0. 

(2) 

Equation (2) signifies that to the left of the dislocation line (x < 0), where the plane 
n = exists, the interaction with the layers n — 1 above and n = — 1 below is the 
same as given by (1). However , to the right of this line (x > 0) the half-plane 
n = is missing and, therefore the layer n — 1 appears to be coupled to the n — —1 
with the same current J c . In our model we ignore any role of the dislocation core, 
because its size is much smaller than the London or Josephson lengths. 

If a vortex resides exactly at and along the dislocation line, the lines of current 
should possess the symmetry with respect to the mirror reflection in the plane z — 0. 
Therefore, no current flows along the half-plane n — 0. It implies that the solution 
for such a vortex in the gauge A z = obeys the conditions: 

2vr 

V$ + — A = , $i = . (3) 

Making use of these conditions, one finds that the potential energy U in (2) should 
be replaced by 

2(l-z/cos%), x<0 

; (4) 

1 — v cos $ s , x > 

Assuming the system is in the 0-SC state {v — 1), the lowest energy solution can be 
achieved for $ s = 0. If, however, v — — 1 (7T-SC), the lowest energy corresponds to 
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the asymptotics 

$ s (-oo) = ±2tt , $ s (+oo) = ±7T (5) 

These account for a half-flux contained in the z = plane. Indeed, employing the 
expression [12] 

? = |£(S.(+oo)-S.(-oo)) . (6) 

for the total Josephson flux ip, one finds that ip = =f(p /2. It is worth noting that 
imposing the condition $ s (+oo) = oo) instead of (5) in order to remove the 
HV makes (4) acquire a contribution proportional to the area either of the half-plane 
n = or the missing part of this plane. Therefore, we conclude that a single edge 
dislocation in 7T-SC binds the HV with an infinite energy. 

To find the energy of the HV in 7T-SC as well as that of the integer vortex in 
0-SC trapped at the dislocation line, we will employ the continuum limit for all the 
layers but the three central ones n = —1,0, 1. It means that the potential energy 
(4) is retained and the summation in (1) is replaced by integration over two regions 
z > +0 and z < — 0. We expand 

s 2 ( 2tt \ 2 

COS ($ n - $ n+ i - Xn,n+l) ~ 1 - -pr [d z <S> + —A z (7) 

2 \ ¥o J 

for the 0-SC and the same for the 7T-SC, with the replacement $ n — > $ n + irn 
made. The continuum approximation turns out to be valid, as discussed in [2,4], 
outside of the junction where the vortex resides, as long as the relation \ a b/s ^> 1 
holds. Therefore, in the following consideration the thickness of the central junction 
made of two infinite layers n — l,n — — 1, with the half-plane n = intervening in 



between, is infinitesimal. As a result, we arrive at the functional (1) written as 

37TV27T/ ,/ 2 <0,2>0 LA^ V </9 7 A^V m Q ) 



87TV27T/ J z <0,z>0 L ^ab ^ ¥0 ' ¥0 

+ ^- f d 2 rdz (curlA) 2 + f d 2 rU 0d , (8) 



k 2 = h -( 2 A 2 j, s 



e Vy? / 

where the conditions (3) and (4), which derives from (3), are taken into account. 

Variation of (8) allows one to obtain the linearized bulk equations (the problem 
is homogeneous in the y-direction and the gauge is A z = 0) as follows: 

2 ^dlA x = ±(d x <S>+ 2 ^A x ) , (9) 
Kb \ <Po J 

and 

2tt 1 

-%,A t = --d,* . (10) 

The solution for the gauge invariant phase = $ + ^r], where the magnetic potential 
A x = d x rj is introduced, can be easily found from Eqs. (9)-(10), if one Fourier 
transforms along the x-axis, yielding 

<M+0) e~ Q " z , z > 

= ; (ii) 

<f) q {-0) e Q « z , z<0 
Q q = \£(l + \y) , 9 (+O) + ? (-O) = . 

Here, <f> q (+0), g (— 0) stand for the Fourier components of the gauge invariant phase 
at the upper and lower edges z = +0, z = —0, respectively, of the central junction. 
Integrating by parts with respect to z in (8) and making use of (9)-(10) one can 
express the free energy in terms of the surface integration only. Taking into account 
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that the magnetic potential rj has no jump at z = and that the phase does have a 
jump, one obtains the line energy (free energy per unit length in the y-direction) as 



e 



i /gov/ L\- _g_ u i 2 

/0 (A r+oo -, ~i 

rfx2(l - z/cos-^) + / dar(l - vcos<p s ) \ , (12) 
-oo 2 Jo J J 



1 

+ 7x* 



<p sq = g (+O) - 9 (-O) . 

It is worth noting that because of the limit \ a b/s 3> 1, the jump of the gauge 
invariant phase (p sq is taken to be equal to the jump $ s (4) of the phase $. In the 
next section we will utilize (12) to calculate the energy of an integer vortex at the 
dislocation line (0-SC case), and show that this energy is less than that of the vortex 
far from the dislocation, implying a binding of the vortex to the dislocation line. 



3 Energy of an integer vortex bound to a disloca- 
tion 

If the penetration lengths X a b were compared to A c , one would have converted (12) 
to the conventional sine-Gordon functional for the JJ [12] making the replacement 
\Jl + \ 2 c q 2 ~ 1, and J2 q q 2 \4>sq\ 2 — * / dx(d x (f) s ) 2 in (12). Consequently, in the 7T-SC 
the HV solution would be obtained. However, the anisotropy in our case will play 
an essential role, implying that (12) cannot be represented in a closed form in the 
physical space of one coordinate x: an infinite number of the gradient terms should 
be retained. To obtain a rigorous upper bound on the vortex energy in the 0-SC we 
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will employ the variational trial function: 

C_e akx + B_e^, 



4> s (x) = < 



C + = 



x < 



(13) 



2tt - C + e 



-kx 



B + e ^ , x > 



1 + a 

a 
1 + a 



2.-(l--)B H 



Here, a, fc, 5± are variational parameters. It can be seen that (13) describes a 
function which is continuous with its first derivative. This form of the trial function 
corresponds to the magnetic field pointing in the y-direction, the total flux being 
<fo (s ee (6)). The term proportional to exp (±x/\ c ) was included to account for the 
fact that at large distances from the vortex the solution of (9-10) decays with such 
exponents along the x-axis. However, it can be shown that all the terms proportional 
to £?'s produce a contribution into the energy as small as s/A a t 1. Thus, in what 
follows, we put B± = in (13). 

Fourier transforming (13) and substituting the result into (12), we arrive at the 
line energy of the vortex residing at the dislocation as 



Bid 
1(a) 



(<Po 


) 2 1 


Utt 




2 




- a7r 






/ a du 


7T 


J 






In {k\ c +ln2-— : 

a 2 — . 

sm-u 2 2 f2(i+ a ) , sinW 



1(a) — — — h O 



s h\f 



ah 



(14) 



u 



+ 



-r 

a Jo 



du ■ 



u 



Minimization of (14) with respect to k gives 



k\ c = /(«)— > 1 

s 



(15) 
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The energy (13) then can be written as 



An) X ab X c 



, X ab _ , n In a 
In — + l + ln2- — 

s a 2 



i + ln/(a) + °(i) 



(16) 



This can be minimized numerically with respect to a, producing 



£ld 



<Po\ 2 1 
An J X ab X c 



ln^ + 1.05 + of-^ 

s \X ab 



(17) 



at a = 0.85, suggesting an asymmetric distribution of currents about the dislocation 
line. 

To find the vortex energy far from the dislocation one can employ the same 
variational function (13), where the symmetric solution (a = 1) should be looked 
for, and the potential energy ought to be replaced by the second line of (4) in a 
whole space as 



U — — J c (1 — v cos 4> s ) , — oo < x < +oo . 
2e 



(18) 



Then, the rigorous upper bound for the integer vortex line energy far from the 
dislocation is 



An) X ab X c 



ln^ + 1.24 + 0f— 

s \X ab 



(19) 



A comparison of (17) and (19) shows that the latter is higher, suggesting that the 
vortex binds to the dislocation. It is worth noting that the exact evaluation of the 
line energy of an interlayer vortex in the ideal layered crystal obtained in [4] gives 



1 \An) X ab . 



abXc 



h ^ + 1.12 + (-£- 

s \X ab 



(20) 



which appears to be higher than the energy (17) of the integer vortex bound to 
the dislocation as well. Comparing (19) and (20), we can say that the varia- 
tional procedure suggested above gives a correct result in the leading logarithmic 
(ln(A a b/s) 1) approximation. At the same time, the second term in brackets 
(19) turns out to be overestimated by 10% in comparison with that in (20). There- 
fore, as soon as the bound vortex line energy (17) contains the same \n(\ ab /c), we 
expect that the binding energy is considerably underestimated. This, however, can 
be improved by introducing additional terms into the variational ansatz (13). 

4 The line energy of the half- vortex in the 7T-SC. 

In the case [8] of the 7T-SC {y = —1) the lowest energy solution obeys the asymptotic 
behavior (5). Therefore, the trial function can be taken in the form 



where the tail ~ exp (±x/A c ) has been omitted, as discussed in Section 3. Calcula- 
tions similar to those performed above for the case v — 1 yield the solution for the 
variational parameters as 




(21) 




a = 0.72 



(22) 



s 



and the half- vortex line energy 




(23) 
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A comparison of this result with (19) shows that in the leading logarithmic approx- 
imation ln(A a6 /s) 3> 1 the line energy of the half-vortex is 1/4 that of the integer 
vortex, as one would expect from estimates [2,4] where the unit flux ip is replaced 
by its half ±(p /2. 

5 Discussion. 

As was emphasized above, an integer vortex turns out to be the same in the both 
states [y = l,u = — 1) of the LD model. It does not allow one to distinguish 
these states by investigating the behavior of such a vortex unless it encounters a 
dislocation. In the 0-SC, where the SC order parameter in the ground state has the 
same sign in all the layers, the edge dislocation binds the vortex with the line energy 
not less than 



This estimate is obtained by a comparison of the exact result (20) obtained in [4] for 
an interlayer vortex energy in an ideal layered crystal with the variational energy 
(17) for the bound integer vortex. 

In the 7T-SC, where the SC order parameters in the ground state have opposite 
signs in adjacent layers [8], the spontaneous half- vortex exists inherently at the dis- 
location line. Therefore, an integer vortex approaching this dislocation will either 
be repelled or attracted to it, depending on the mutual orientation of the two vor- 
tices. If their moments are opposite, the integer vortex will recombine with the 




(24) 
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half-vortex. As a result of such a process, the integer vortex disappears, and the 
half-vortex switches orientation. It means that the energy of the former will be 
released in some way. 

Given these two different types of interaction of free vortices with a linear edge 
dislocation, the two possible states of the LD model can be distinguished by inves- 
tigating the dynamics of vortices in a single HTS crystal with a density of interlayer 
edge dislocations less than the reciprocal area (l/A a feA c ) occupied by one vortex. 
Under this condition the pairs of dislocations and vortices bound to them can be 
considered independently of each other, and our results will apply. Thus, in the case 
of the 0-SC these dislocations should play a role of weak pinning centers. Therefore, 
after imposing and then removing an external magnetic field along the layers some 
residual magnetisation might be observed, because of the vortices trapped by the 
dislocations. Effects of thermal fluctuations and weak vortex-vortex interactions will 
cause this magnetisation to decay slowly with the typical time determined by the 
line binding energy (24). 

As it can be shown, in the 7T-SC the ground state of an array of half-vortices 
attached to dislocations in a single crystal will be antiferromagnetic. Hence, no 
macroscopic magnetisation is expected unless an external field along the layers is 
imposed. Introducing external vortices, with their density being twice less than 
that of the dislocations, will result in them recombining with that part of the half- 
vortices which have the orientation opposite to the external field. As a result, a 
residual magnetisation of the sample will occur. This magnetisation, which was 

12 



originally due to the integer vortices, is now produced by the half-vortices ordered 
ferromagnetically. With the external field removed, such a configuration is not the 
lowest in energy, and these half-vortices should re-orient themselves to return to 
their former antiferromagnetic ordering. However, to do so the half-vortex should 
emit an integer vortex which, then, is to be expelled from the sample. The energy 
required for such a process is that of the integer vortex , and, as comparison of 
(19), (20) with (24) shows, it is bigger than the binding energy (24) in the 0-SC by 
the factor of ln(X ab /s) ^> 1. Therefore, in the 7r-SC the residual magnetisation, 
which is controlled by the dislocation density, should decay with a typical time 
constant much bigger than in the 0-SC. 

The other implication of the 7T-SC was discussed in [11]. If a SC grain is in the n- 
SC and contains an odd number of dislocations, it should carry a net magnetisation 
±(/? /2 because of the uncompensated half-flux. An array of such grains coupled 
by means of the magnetic field only would exhibit paramagnetic response. In this 
respect the paramagnetism of the HTS granular materials [13] might be considered 
as an indication of the existence of the 7T-SC. The recent observation [14] of the 
half-flux trapped inside a HTS ring could also indicate that the grains in this ring 
are in the 7T-SC, and that they are joined in such a way that the Cu-0 layers form 
a screw-like structure, while traveling around the ring. In this case such a ring 
characterized by the structural chirality should develop a spontaneous half-flux. 

In conclusion, we employed a variational approach to show that a single interlayer 
edge dislocation binds an integer interlayer (coreless) vortex, within the framework 
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of the LD model. The current distribution around such a bound vortex turns out 
to be asymmetric. Regarding the suggestion [8] of the 7r-type interlayer Josephson 
coupling, we have shown that half-vortex is attached to the dislocation line, and 
calculated the half-vortex line energy. The physical consequences of such a ground 
state were discussed. 

We are grateful to John R. Kirtley for useful discussions of possible methods for 
experimental detection of half-vortices. 
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Figure Caption 



Fig.l Schematics of the edge dislocation made of half of the conducting plane 
n = (z = 0,x < 0, — oo < y < +00 ) which is inserted between the two layers 
n — ±1. The axis z is perpendicular to the layers far from the dislocation core, 
which location is indicated by the vertical bar. 



16 



